Sr2(Ba2)V04 under pressure — an orbital switch and potential superconductor 
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We study Sr2(Ba2)V04 under high pressure by means of the local density approximation + 
dynamical mean field theory method. While Sr2V04 is a 1/6-filling three-band system at ambient 
pressure with a small level splitting between the dxy- and dy^/ ^^-hands, we show that an orbital 
polarization occurs under uniaxial pressure, resulting in dramatic changes of the magnetic, optical, 
and transport properties. When pressure is applied in the c-direction, a analog of d^ cuprates is 
realized, making Sr2(Ba2)V04 a possible candidate for a superconductor. Experimentally, this 
uniaxial pressure can be realized by growing Ba2V04 on a substrate with lattice constant 4.1-4.2 A. 

PACS numbers: 71.27.-|-a, 71.30.-|-h 



Since the discovery of high temperature superconduc- 
tivity in cuprates [H, strongly correlated electron systems 
(SCES) and their intriguing magnetic, dielectric, optical 
and transport properties have been at the center of solid 
state research. Hence, a quantitative reliable calculation 
of correlation effects from first-principles is one of the 
most important challenges. This is particular difficult 
since the standard local density approximation (LDA) 
in the framework of density functional theory[2|] fails if 
electronic correlations are strong. Recently however, a 
variety of attempts which go beyond LDA have been un- 
dertaken, and many successes have been achieved jsj. 

The next step in this direction is the (theoretical) ma- 
terials design of SCES with specific properties and the 
simulations of SCES under extreme conditions. As a 
touchstone for such attempts, we study the electronic 
structure of Sr2V04 and Ba2V04 under high pressure by 
means of the LDA+DMFT (dynamical mean field the- 
ory) method^, one of the most widely used approaches 
for realistic calculations of SCES ^] . 

The reason why we focus on Sr2V04 is twofold. 
First, Sr2V04 is a layered perovskite, as cuprates and 
ruthenatesQ which show unconventional superconduc- 
tivity. Second, the challenge to synthesize single- 
crystalline Sr2V04 has been overcome quite recently: 
One of the authors and his coworkers [tJ employed epi- 
taxial growth techniques for growing a thin Sr2V04 film. 
Hence, a detailed investigation of the electronic structure 
becomes now possible. 

As for the first point, in fact, Sr2V04 attracted atten- 
tion because it is a "dual" material of La2Cu04. Namely, 
the former has one 3d electron per V site {d^ system), 
while the latter has nine 3d electrons per Cu site {dP sys- 
tem). However, as was already pointed out by Pickett et 
al. in 1989 Q, there is a big difference between these ox- 
ides. While La2Cu04 is a 1/2-fining single-band {dx^^y2) 
system, Sr2V04 is a 1/6-filling three-band {dxy/yz/zx) 
system, since the level splitting of t2g in the latter mate- 



rial is much smaller than that of Cg in the former. There- 
fore, the theoretical idea of unconventional superconduc- 
tivity in Sr2V04 has been dismissed. 

Also experimentally, Sr2V04 and La2Cu04 behave in- 
deed differently. Magnetic properties and transport prop- 
erties were rneasured for a polycrystal of Sr2V04 in the 
early 1990'si9|, indicating that Sr2V04 is an antifer- 
romagnetic insulator (semiconductor) with a low Neel 
temperature T/v ~ 45K. But in contrast to La2Cu04, 
a small ferromagnetic moment was also observed. In 
recent measurement of the optical conductivity for a 
single-crystalline thin film, a small gap structure, i.e., 
a peak around 1 eV and a shoulder around 0.5 eV was 
observed [?!. 

Theoretically, several first-principles calculations be- 
yond LDA were performed, and it was confirmed that the 
relation between Sr2V04 and La2Cu04 is not dual|l3, 
11]. Especially, Imai et al. carried out an LDA-I-PIRG 
(path integral renormalization group) calculation and 
found a nontrivial orbital-stripe order [lO|. These orbital 
degrees of freedom are irrelevant in cuprates. 

The motivation of the present study is based on the fol- 
lowing idea: While Sr2V04 and La2Cu04 are certainly 
not dual at ambient pressure, we hope to change this 
by changing the atomic configration of the octahedron 
around the V ion. Here we consider to apply high uni- 
axial pressure, synthesize films on substrates with appro- 
priate lattice constant, and introduce chemical pressure 
(substitute Sr by Ba). If the level splitting between the 
dxy and dy^/zx orbitals becomes larger, it will lead to an 
orbital polarization, so that the system might actually 
become a d^ analog of d^ cuprates. In this Letter, we 
examine this idea by means of LDA-I-DMFT. 

GGA optimization of crystal structure. Let us now 
turn to our actual calculations. First, we perform a 
GGA (generalized gradient approximation) calculation 
with plane-wave basis, ernploying the Tokyo Ab-initio 
Program Package (TAPP)[lJ]. Note that a plane -wave 
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1: (a) Atomic configuration of Sr2V04 togetfier witfi tfie 
NMTO Wannier function of tfie V d^^y-orbital. The position of 
Sr and the apical oxygen are optimized by a plane-wave GGA 
calculation, (b) GGA-optimized distance between O and V 
along the c direction (open squares) and in the ab-plane (open 
circles) as a function of the c axis elongation c/co- 
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basis set has advantages for optimizing atomic configura- 
tions under high pressure which are unknown experimen- 
tally. We adopt the exchange-correlation functional in- 



troduced by Perdew et al.jl 
tials in a separable form[l 




and ultra-soft pseudopoten- 
ISl. The wave functions are 



expanded up to a cut-off energy of 36.0 Ry, and 8x8x8 
and 12 X 12 X 12 fc-point grids are used. We assume that 
the system has the same lA/mmm symmetry as at at- 
mospheric pressure [l6| . so that there are only two free 
parameters, i.e., the position of the apical oxygen and 
Sr. In Fig.[TJa), we show the atomic configuration of the 
octahedron around the V ion. 

We calculate the total energy as a function of the lat- 
tice constant a, fixing the ratio c/a to its experimental 
value (=3.28) [9j. We find that the energy minimum is at 
a=3.89A(not shown), in excellent agreement with exper- 
iment (a=3.84A). When the lattice constant a is fixed to 
this optimized value, the Sr-V and 0-V distances along 
the c-direction become 4.46A and 2.0lA, respectively, 
consistent with the experimental values 4.46A and 1.98A. 

Next, we change the lattice constants to simulate the 
effect of pressure. Considering to apply high uniax- 
ial pressure or, more realistically, to synthesize Sr2V04 
films on substrates with appropriate lattice constant, we 
change the lattice constant c up to ±10% of the experi- 
mental value Co =12.76A. 

In Fig. [TJb) , we plot the GGA-optimized distance be- 
tween O and V along the c-direction {dc) and within the 
a6-plane (dab) as a function of c/cq. Here we fix the vol- 
ume of the unit cell (V) rather than a, since the total 
energy at fixed V is always lower than that of fixed a. 
While the V-0 distance dc is longer than dab at ambient 
pressure, dc becomes shorter than dab for c/cq < 0.95. 

In fact, the ratio dc/dab determines the splitting of 
the three t2g orbitals. If we press uniaxially along the 
c-direction, the negatively charged oxygen ions move to- 
wards the vanadium site. Hence the energy of the dyz/zx~ 
orbitals, which point along the c-direction, is enhanced. 
At ambient pressure the level splitting is small. But given 



FIG. 2: Density of states of the t2g band for (a) atmospheric- 
pressure, (b) uniaxial pressure in the c direction (c/co = 0.9), 
and (c) uniaxial pressure in the ab plane (c/co ~ 1.1). Solid 
(dashed) line is for the d^jz-orbital (d^^/^^-orbitals). Inset: 
LMTO band structure for atmospheric-pressure. 



the fact that dc and dab change considerably with c/cq 
in Fig. [TJb) , we may expect that we can control the level 
splitting, and, consequently, the orbital occupation by 
applying pressure. 

LDA+DMFT calculation. To examine this idea 
in the presence of electron correlations, we perform 
LDA+DMFT calculations for the atomic configurations 
obtained above. To this end, we first carry out LDA 
band structure calculations with the LMTO (linearized 
muffin tin orbital) basis 17[ ■ In the inset of Fig. [2l we 
show the obtained band structure for ambient pressure. 
Almost the same band structure is obtained by the GGA 
calculation with plane- wave basis. 

Then, we extract the three t2g bands by the NMTO 
(7V-th order muffin tin orbital) downfolding[18|, using the 
generated LDA potential. As a typical example of the 



xy 



resulting NMTO Wannier functions, we show the d. 
orbital of V in Fig.[T](a) for ambient pressure. The density 
of states (DOS) of the t2g band for c/cq = 0.9, 1.0 and 
1.1 are shown in Fig. H The NMTO band width of the 
dajy-band is 2.46 eV at ambient pressure (c = cq). This 
value is consistent with that of the GGA calculation with 
plane-wave basis and that of Pickett ct al. who used the 
full-potential linearized augmented plane wave (LAPW) 
basis 0. Imai et aL[l3] reported a smaller band width of 
-2.0 eV0. 

For the case of c/cq^O. 9(1.1), Fig. [5] shows that 



the 



center of gravity of the dyz/zx 
higher (lower) than that of the dxy-h&nd 



bands is clearly 
Indeed, the 



crystal field splitting between the dxy- and d^^/yz- 
orbitals in the NMTO Hamiltonian is -382(-h434) meV 
for c/co=0. 9(1.1). Concerning the electron occupation of 
the dxy- and rfjj^/jjx-orbitals, 90%(20%) of the d electrons 
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FIG. 3: LDA+DMFT spectral function for atmospheric- 
pressure, where dashed, dotted, and sohd hne are for dxy, 
dxz/yz, and the total Sr2V04 spectrum, respectively. For 
U' = 2.8, two solutions [metallic (b) and insulating (c)] co- 
exist. 



are accommodated in the d^jy-band for c/co=0. 9(1.1), in 
contrast to ambient pressure where all three bands 
are similarly occupied. 

Next, we perform DMFT calculations for the three low- 
energy t2g bands, studying whether electronic correla- 
tions result in a full orbital polarization for c/co=0.9 and 
1.1. To this end, the DMFT effective impurity model is 
solved by the standard Hirsch-Fye quantum Monte Carlo 
(QMC) method [20|. where the temperature is 0.1 eV with 
100 imaginary time slices and the number of QMC sam- 
ple is ~ 2 X 10^. We employ the relation U = U' + 2J 
where U, U', J are the intra-orbital Coulomb interaction, 
the inter-orbital Coulomb interaction and the Hund cou- 
pling, respectively. 

We first calculate the spectral function for ambient 
pressure with various interaction parameters C/', fixing 
J = 0.7 eV which is a reasonable value for V ions. In Fig. 
131 we show the spectral function for U' = 2.5 cV (a) and 
2.8 eV(b,c); the Mott-Hubbard transition occurs around 
U' = 2.5 ~ 2.8 eV. Since the gap in (c) is in accord with the 
main optical peak in experiment we expect f7'^2.8cV 
for Sr2V04. The coexistence of this insulating solution 
(c) with a metallic one (b) indicates insulating Sr2V04 
to be close to a Mott-Hubbard transition. 

Note that Imai et al. [l3| estimated smaller Coulomb 
interactions: U' ~ 1.3 and J ~ 0.65 eV by the combi- 
nation of the constrained LDA and the GW method [13] ■ 
However, as mentioned above, the band width of the ^29- 
orbitals of Ref. [l3| is 20-25% smaller than that of the 
present study. If we normalize the interaction parame- 
ters by the band width, the difference is not so big. On 
the other hand, Sekiyama et al. employed U' = 3.55 and 



J = 1.0eV in their LDA-f-DMFT calculation for SrVOa, 
reproducing the photoemission spectrum of SrV0^[2l|. 
These values are not so far from ours. 

Let us now turn to the LDA-I-DMFT results for high 
pressure. In Fig. HJa), we plot the spectral function for 
c/cq = 1.1, showing a metallic peak at the Fermi level. 
In contrast to ambient pressure, there is no coexisting 
insulating solution at c/cq = 1.1, i.e., applying pressure 
in the ab plane makes Sr2V04 metallic. An important 
point is that the d^^y-orbital is almost empty for c/cq = 
1.1: orbital polarization occurs. The system becomes a 
quarter-filled 2-band Hubbard model which is well known 
to have a ferromagnetic ground state 2^, 2^- Therefore, 
we expect ferromagnetic spin fluctuation to be dominant 
at low T if pressure is applied in the ab plane. 
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FIG. 4: Same as Fig.[3]but for (a) uniaxial pressure in the ab 
plane (c/co = l.l) with U' =2.8 and (b) uniaxial pressure in 
the c direction (c/co = 0.9) with U' = 2.5. 

On the other hand. Fig. |DJb) shows the result for uni- 
axial pressure along the c direction (c/cq = 0.9). We see 
that the spectrum is now clearly insulating, even for the 
smaller value of U' = 2.5 for which we have a metal in 
Fig. [3ja). This is a surprising result: an inverse Mott- 
Hubbard transition. (Usually, applying pressure makes 
an insulator metallic, not a metal insulating as from Fig. 
[31[a) to Fig. HKb).) What is the reason for this unusual 
behavior? From Fig. [21 we see that the LDA band width 
does not change strongly from 2.26 eV at c/cq = 1 to 2.46 
eV at c/cq = 0.9. This small change of band width alone 
would indeed indicate more metallic behavior-as usual. 
But more important is that uniaxial pressure changes the 
crystal-field splitting: The two dy^/^^'Orbitals become 
unoccupied, and then the large intra-orbital repulsion U 
makes the remaining (single) c?a;y-orbital Mott-insulating. 
This way, a analog of cP cuprates is realized. Concern- 
ing the magnetic properties, we expect an antiferromag- 
netic instability since the system becomes a half-filled 
single-band model. 

502^04. While it is an interesting possibility to 
change the electronic properties by controlling the lat- 
tice constant, it might be difficult to change its value up 
to ± 10%. For exapmle, in order to grow Sr2V04 thin 
films with c/co = 0.9, we need considerably large lattice 
mismatch between a substrate and bulk Sr2V04, while 
an excessive mismatch often results in lattice relaxation. 
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FIG. 5: Density of states of the t2g band of Ba2V04 for 
(a) atmospheric-pressure and (b) uniaxial pressure in the c 
direction (c/co = 0.95). Figs, (c) and (d) are LDA-I-DMFT 
spectral function for (a) and (b) respectively, with U' = 2.5. 



Thus, lastly, we consider the possibility of chemical 
pressure by substituting Sr by Ba. Since the ion radius 
of Ba is larger than that of Sr, the crystal is expected 
to be expanded. The important point here is that the 
V-0 distance in the ab plane and that along the c axis 
will increase differently. Namely, while the increase of a 
directly affects the V-0 distance in the ab plane (the lat- 
ter is exactly half of the former), c and the V-0 distance 
along the c axis are independent parameters. 

First, we perform GGA calculation with structure op- 
timization for Ba2V04 by changing a and c. The energy 
minimum is at a=4.04 A and c/a — 3.36, the V-0 dis- 
tances are 2.02 A in the ab plane (compared to 1.92 A 
for Sr2V04) and 2.01 A along the c direction (2.01 A for 
Sr2V04). This means that the situation for Ba2V04 is 
similar to that of Sr2V04 with c/cq — 0.95, i.e., if uniax- 
ial pressure is already applied. 

Next, we consider the compression c/cq — 0.95 for 
Ba2V04. The orbital polarization becomes larger, i.e., 
the electron density is 0.73 for the dxy band and 0.14 
for the dyz/zx band. Thus it is interesting to proceed 
with DMFT calculations. In Fig. [5l we plot the result- 
ing LDA-I-DMFT spectral function for ambient pressure 
and uniaxial pressure in the c direction (c/cq = 0.95) with 
U' = 2.5, along with the density of states by LDA. We 
can see that a large orbital polarization is realized even 
for c/cq = 0.95. Since the relative hopping and inter- 
action parameters for the d^y band are very similar to 
the cuprates, a d^ analog of d^ cuprates can be realized. 
Experimentally, the uniaxial pressure is best realized by 
growing Ba2V04 on a substrate with a larger in-plane 
lattice constant 4.1-4.2 A. 

We thank J. Akimitsu, O. K. Andersen, H. Aoki, H.-U. 
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